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ABSTRACT 

' We present a new class of exact nonsingular solutions for the Maxwell equations in vacuum, which describe the 
I electromagnetic field of the counterpropagating focused laser beams and the subperiod focused laser pulse. These 
' solutions are derived by the use of a modification of the " complex source method" , investigated and visualized. 
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1. INTRODUCTION 



00 ' According to predictions of QED, vacuum can be modified by strong electromagnetic fields if the strength 
becomes greater than Fq = m?c^ /eh = 1.32 ■ lO^^V / cm and, in particular, it should behave as a nonlinear 
media with respect to transmission of such fields. Different aspects of this phenomenon (including, e.g., vacuum 
C/3 , polarization, elastic photon-photon scattering, vacuum birefringence, photon splitting, pair production, etc.) had 
. ^ ' been investigated theoretically very actively, especially after the beginning of sixties. However, at that time the 
] researchers restricted themselves by considerations of the simplest configurations of the electromagnetic field, 
Q namely, a constant homogeneous field and a plane wave field. It was motivated, partially, by simplicity (the Dirac 
^ , equation in these backgrounds admits exact solutions). But also, it was commonly doubted that the fields with 
Q • the strength of the order of Fq will be really created in a laboratory in viewable times. However, it was declared 
' ^ ' recently that this experimental task could be handled on the basis of the modern level of technology^"^ and 
moreover, resolution was announced to happen already in the coming decade.^ The high-intensity field, which 
' can be obtained by implementation of the proposals cited above, will be most probably realized in ultrashort 
(maybe, subperiod) tightly focused laser pulses of optical range, though alternative projects based on XFEL are 
also under consideration. Therefore, in order to plan experiments in a near future, one should, firstly, examine 
T-H ■ whether previous calculations of quantum effects in a strong field can be applied to more realistic fields, and 
^ ' secondly, look for those new effects, which can arise in realistic fields but vanish identically in the homogeneous 
' or the plane wave fields. 

' In view of this, study of nonlinear QED effects has being risen essentially. Let us mention here just a 

few of these recent studies, which are most closely related to consideration of the realistic structure of the 
I/"} ' electromagnetic fields expected in future experiments (for a review of other popular experimental proposals 
' see*). In the papers^"^ scattering of relativistic electrons on a focused laser beams and pulses was considered. 
In particular, in^'^ the asymmetry of the scattered electrons observed on the experiments was explained. The 
papers*'^ dealt with pair creation from vacuum by a super intense (E Fq) focused laser pulse or by two colliding 
pulses. In Refs.^"'^^ a new effect of odd harmonics generation by a tightly focused laser beam in vacuum was 
investigated. Particularly, it was shown in Refs.,^' ^"^^ -'^^ that the vacuum quantum effects induced by a tightly 
; ^ focused laser beam or pulse can become observable when the peak strength of the field will reach the values of 
magnitude of about one order less than the characteristic value Fq. Another message of these papers was that 
typically these effects are crucially sensitive to the parameters and the particular type of the model of the field. 

However, consideration of the studies mentioned above was based on the approximate model of the focused 
laser field, which was suggested in Ref.^ (consideration of Ref. was based on a similar but quite different model 
which was suggested independently by K. McDonald though not published). This model can be considered as a 
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generalization of the well-known paraxial approximation for a Gaussian beam to the case of a transverse vector 
field (note, that in Ref.^ , among the other issues, the classification of polarization types of the focused beam 
was discussed exhaustively). To take into account finiteness of temporal duration of the pulse, one, roughly 
speaking, multiplies the carrier field by the slowly varying temporal envelop. It is obvious that this approach is 
valid only for sufficiently long and weakly focused pulses. In order to check the validity of these approximations, 
it is desirable to have exact focused beam and pulse-like solutions of the Maxwell equations in vacuum. 

Perhaps, the most elegant attempt to construct exact analytical beam or pulse-like field models was made 
in Refs.^^"^^ for a scalar field and was recently directly generalized to the case of transverse vector fields in 
Refs.^^'"'^^ The main idea of this technique (now called in the literature "the complex source method") is to 
generate the focused fields in vacuum by the pointlike sources whose position is shifted into the complex domain 
(this is mathematically equivalent to analytical extension of a retarded potential). A first exact analytical model 
of the focused electromagnetic pulse obtained by this method was manifested in Ref.^* . Being compared with 
the approximations described above, this model revealed two new qualitative features. First, the time average of 
the electric field at any fixed spatial point vanishes identically, which is very reasonable. And second, it possesses 
some temporal analog of the Gouy phase which provides chirping at the wings of the pulse. This last effect was 
called in " the self-induced blueshift" . 

In this paper we show that the usual considerations of this complex source method suffer from two disad- 
vantages. On the one hand, the focused field models presented in the literature thus far^^"^* possess algebraic 
singularity, which can be attributed to some real sources (unlike the fictional pointlike "complex sources") lo- 
cated on a ring in the focal plane. Therefore, in a sense, these solutions simulate propagation of the focused 
field through a circular aperture bounded by the branch cut in the focal plane, rather than the field focused in 
vacuum. Only in the special limiting case of weak focusing (when the paraxial approximation works well) this 
singularity escapes away to infinity and does not affect the field in the focal spot. On the other hand, as we 
show below, the complex source method provides analytical solutions only in the form of two colliding beams or 
pulses (in the limit of weak focusing one of them disappears). We also sketch the method for separating a single 
pulse solution and discuss the possibility of its representation in a simple analytical form. The consideration 
below is arranged as follows. In the Sec. [5] we remind how a certain class of electromagnetic field configurations 
can be parameterized by an axially symmetric scalar field configuration satisfying the usual wave equation. This 
transformation is very close (in fact, coincides) to the consideration of Ref.^ However, in order to increase the 
clarity of presentation we use different notations and reasoning. Existence of this mapping between the scalar 
focused fields and a wide class of transverse vector focused fields with clear polarization type (say, circularly 
polarized e- or h-fields) prevents us from consideration of the more complicated vector version of the " complex 
source method" Therefore, in Sec. [3] we restrict ourselves to the more clear scalar version of the complex 

source method and demonstrate its drawbacks. The modifications which resolve these problems and the non- 
singular beam- and pulse-like solutions of the Maxwell equations are presented in Sec. [H Our conclusions are 
collected in Sec. [H 

2. THE NAROZHNY-FOFANOV MODEL 

Consider a laser beam or pulse focused along the z-axis. In the three-dimensional transverse gauge, the electro- 
magnetic field is expressed through the vector potential A, which satisfies equationsQ 

A-AA = 0, V-A = 0. (1) 

Let p, (j) and z be the cylindrical coordinates. Since the pulse is confined in radial direction, the electric and 
magnetic fields are generally not orthogonal to the direction of propagation of the focused beam or pulse. 
However, as usually, an arbitrary field can be represented as a superposition of the e-field (E orthogonal to 
z) and the h-field (H orthogonal to z). Since an h-field can be always obtained from the e-field by a duality 
transformation E ^ H, H ^ — E, it is enough to consider, say, the case of the e-field. Hence, in the following 
we assume Az = 0. In the simplest case of circular polarized field, the pulse consists of photons with jz = +1, 

*We use such units that c = ft = 1. 
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where j is the total angular momentum of the photon (assembled from the orbital momentum 1 and the spin s). 
This means, that 

ISA . , . 

-— +zzxA = A, 2 

I a<p 

where z is the unit vector in z-direction (compare to discussion of spherical waves in Refs.^^'^°). The general 
solution to ^ can be written in the form 

A(p,0,z,t) = e+Fi(p,z,t)-e_i^2(/o,^,t)e2'^, (3) 

where e± = x±iy, and Fi^2 ~ the arbitrary functions (evidently, the latter equation provides decomposition of a 
state with jz = Iz + Sz — +1 onto the states with Iz = 0, Sz = +1 and with Iz = +2, Sz — —1)- By substitution 
of the anzatz ([3]) into the equations ^ one obtains the following equations for the auxiliary functions ^1,2, 

AFi = i {pFip)p + Fiz. = Fi , (4) 
A - ^F^ = F2, (5) 

Fip = F2p + -^F2. (6) 

The electric and magnetic fields can be expressed in terms of functions F12 by 

E = -A = - (e+Fu - e^F2te'''*') , (7) 
H = V X A = (e+Fi^ + e^F2ze^"'' - 2 zFipe*"^) . (8) 

The equations ^ and ^ are compatible and F2 can be expressed through Fi by 

p 

F2 = Fi~^ j Fipdp. (9) 


However, it is more suitable to express both functions Fi^2 in terms of a unique auxiliary function U{r,z,t) 
through 

Fi,2 = Upp ± -IJp, (10) 
then the equation ([H]) is satisfied identically. The new function U obeys just the scalar wave equation 

C/-Af7 = 0, (11) 
and it is easy to check that the vector potential can be expressed in the following compact form 

A = 2[e+ ■ (z X V)](z X V[/). (12) 
For example, by choosing an approximate monochromatic solution (i.e., applying a paraxial approximation) 

for the wave equation pT|) . where Eq is the field strength amplitude, ^ — p/R^ x = z/L, R- the focal radius 
(i? ^ A = 2tt/lo), L — ojR^- the diffractive length, and (p — u){t — z)- the fast phase, one comes to the particular 
focused beam field model, which was previously proposed by Narozhny and Fofanov and successfully applied to 
the problems of electron scattering,^' ^ pair creation^' ^ and harmonics generation^"' by a focused laser field. 

In this way it follows that in order to construct a realistic, vector model of the field in a focused electromagnetic 
beam or pulse, it is sufficient to consider a cylindrically symmetric scalar model first, and then apply the mapping 
(fT^ . In the next section, let us consider the original complex source method for derivation of scalar focused 
beam (pulse)-like solutions of the Eq. ([TT|l . 
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3. THE COMPLEX SOURCE METHOD 

The idea of the complex source method looks as follows. First, consider a scalar wave equation with an arbitrary 
variable pointlike source q{t) located at the origin, 

U - AU = 4:nq{t)S{r). (14) 

Obviously, the created field is expressed in terms of the usual retarded potential, 

+00 

f/(r,t)= / D^'-''\v,t~t')q{t')dt' ^ (15) 



where (r, i) is a retarded Green function for the wave equation. 

Now, let us shift the location of the source into the complex domain of the z axis, r — > r — ihi where 6 > 0. 
It is also suitable to accompany this by a temporal shift t t ~ ib. Then, the analytical continuation of the 
retarded potential (fT^]) reads 

U{r,t)^ ''^'~'^-'^\ n'=p^ + {z~^hf. (16) 

In addition, let us choose the branch of the square root in the definition of TZ so that Im(7?.) < 0, then in terms 
of the standard root we have TZ = {z — ib)^Jl + / {z — ib)"^. Then, according to the arguments presented in 
Refs.^^"^^ , the real part of the field ([TC]) describes a scalar field focused along the z axis and, in particular, 
satisfies the source free wave equation (jlip. The viewable meaning of this complex shift transformation can be 
better understood by consideration in the k-space. The equation ()15p can be represented in the form 

47r 

^•'^ ^ "(c. + zO)2-P'^"' (^^^ 
and the shift transformation is reduced to Uki^ — > cxp[b{kz — Lu)]Uki^. On the mass shell, we have w = fc so that 

Uk^ ^ exp[-&cj(l - cos6')][/kc^ « exp(-6cj6'V2)[/kc., b > A (18) 

where 9 denotes the angle between k and the z axis. 

Thus, the complex shift suppresses those wave vectors in Fourier decomposition of the radiated field which 
are highly inclined with respect to z, as it is indeed expected for a focused field. For example, if one chose 
q{t) = goe"*"* then, in the limit 6 > A, p < (6A)^/^ < (6^A)^/'^ by applying the expansion TZ « z-ib+p'^ /2{z-ib), 
one can reduce ([T6l) to the form (fT3|) with L — 2b and Eo = —Suj'^qo/L'^. On merely these grounds it was proposed 
in the literature to introduce an additional envelop factor in q{t) and to consider the equation (I16p as an exact 
focused pulse-like solution of the wave equation with the shift parameter b identified with the half of the diffractive 
length. 

Nevertheless, the above arguments are not convincing enough for the following two reasons. First of all, an 
attempt to consider the Eq. (|16p as a solution for the homogeneous wave equation is incorrect and is based on the 
illusion that the source term in the RHS of the Eq. (fTi|) after the shift, being proportional to 5{r — ibz), should 
vanish identically when the vector r possess real valued components. In fact, a (5-function of a complex argument 
behaves in a more complicated way than its real argument analog. Presence of the sources follows already from 
a simple observation that the transformation (fT5)) does not replace the pole singularities in the expression ([T7]) 
at the mass shell by 5{u!'^ — fc^). The location of the sources coincides to the ring singularity p = b, z = oi 
the expression (11611 1. This ring is picked out physically by zero delay time for radiation emitted by the complex 
source and crossing the physical hypersurface Im(z) = 0. Moreover, when approaching this singularity along the 
focal plane z = we have TZ oc \/\p — b\, so that U oc \p — b\^^/^ . This growth is faster than log[l/|p — b\], which 



^Strictly speaking, due to untypical choice of the branch of a square root for TZ, the sources are present in the exterior 
of this ring at the focal plane, which is the branch cut, as well. 
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could be expected if the linear density of sources was finite at the ring. Since the focal radius R ^ (6A)^/^, the 
singularity becomes safe due to its escape sufficiently far away from the focal region which contains the most 
part of the field energy only in the weak focusing limit b ^ X. However, the paraxial approximation works well 
exactly in this limit, so that knowledge of the exact solution becomes unnecessary in this case. 

Another drawback of the Eq. ^TE\i can be seen from Eq. ^TE\\ . Indeed, although due to a complex shift the 
wave vectors with 6 3> {buj)~^^^ become suppressed relatively to those with 6 < (fow)"^/^, nevertheless they 
remain be present. The factor exp[5(fc — fc^)] just mimics but can not replace the true factor 9{kz) which is 
expected to enter the Fourier decomposition of a single beam or pulse focused along z. Therefore, the Eq. (fTH|l 
actually describes two counterpropagating focused fields of non-commensurable amplitudes. Only in the limit 
6 ^ A the amplitude of a pulse propagating in the direction opposite to z becomes exponentially suppressed. 

4. THE NONSINGULAR BEAM AND PULSE-LIKE SOLUTIONS 

The most obvious and direct way to correct the Eq. and obtain the true, nonsingular solution for Eq. PTI) 
is to take a difference of the retarded and advanced solutions, i.e. to start from 



+ 00 



t/(r,i) 



(r, t-t')- D'-'"^"^ (r, t - t')] q{t') dt' 



q{t - r) - q{t + r) 



(19) 



instead of (fT5|) . Being a difference between two solutions of the same inhomogeneous equation, the field (fT9|) 
undoubtedly represents a solution for the corresponding homogeneous, i.e., the wave equation. In particular, 
pop is non-singular at the origin (if q is differentiable, of course) and its Fourier transform 



-47r 



87r^isgn(w)(5(tJ^ - k'^)qu. 



vanishes everywhere except on the mass shell. 

Applying the complex shift transformation of the previous section to solution p9p . one obtains 

q{t-ib-n) -q{t-ib + n) 



C/(r,t) 



n 



n=^p^ + iz- ^fe)2, 



(20) 



(since (|19p is even with respect to r, an arbitrary branch of the square root defining TZ can be chosen in this 



In order to construct a monochromatic focused beam-like solution, let us choose q{t) — q^e where qo is 
a normalization constant. Then we have 



U{r, t) — 2iqo exp[— tj^ — iLot] 



(21) 



Particularly, in paraxial approximation \, p < (foA)^/^ ^ (6^A)^/^], by expanding 71^:1 z — ib + /'2,{z — ib), 
we obtain 



C/(r,i) 



b{l + izlb) 



exp 



iuj[z — t) 



ujp 



2b{l + iz/b) 



exp 



-2Lub - iLu{z + t) + 



ujp 



2b{l + iz/b) 



The second term in the parentheses (which had arisen due to a modification proposed here) is exponentially 
suppressed since uj^^, p <^b. Therefore, by comparison with Eq. (|13p . we can identify the parameter b with a 
half of the diffractivc length L and go with —EqL'^/'&uP'. In the opposite (formal) case of ultra tight focusing 
6 <C A, our solution (j2ip reduces to a superposition of smooth convergent and divergent spherical waves, unlike 
()16p . which becomes a convergent spherical wave in the half space z < and a divergent wave in the half space 
z>0. 
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The auxiliary functions F12 introduced in Sec. [21 in the case under consideration read 

Fi = ^^<loeM^^b-iujt] ^ ^2^2 _ 3^2] _ sin(c.7^)] - iu^p^U' sin(c.7^) } , (22) 

F2 = _ 2»goexphu;&-zu;t] |3^^^^g(^^) ^ ^^2^2 _ 3) 31^(^7^)} , (23) 



and the components of the vector potential and the electric and magnetic fields can be obtained by plugging the 
expressions ((^ . ((^ into the formulas Eq. ([7]) and ^ and taking the required derivatives. 

As well as the ordinary complex source method, up to now our modification always generates pairs of coun- 
terpropagating focused fields. In order to separate out a single focused beam, one can multiply the Fourier 
transform of the field by the factor 9{±kz), which is equivalent to a transformation U U±, 

+00 
—00 

The integral in (|24p can be evaluated by enclosing the integration contour and computation of residues if U 
vanishes at infinite z, e.g., if it is given by a rational function. In this particular case, the transformation 
projects U onto a space of functions which are analytical in the upper half of the complex plane z. 

Let us discuss in brief application of this transformation to the beam-like solution (PTjl . For this purpose, let 
us first represent it in a more viewable explicit form 



exp 
U (X 



iLu(z ^ ib)^l + —exp —iLu{z — ib)^^JT+ 



(25) 



For definiteness, in what follows consider extraction of U+. After substitution into the formula ([M]) one obtains 
two integrals over z', each one corresponding to one's exponent on the RHS of Eq. (PS]) . The first integral can 
be enclosed in the upper half plane, while the second integral - in the lower one. Furthermore, in the case 
p < b the first integrand in the upper half plane possesses a pole at z' = z + iO and the two branch points at 
z' = i{b±p), while the second one is analytical in the lower half plane. Consequently, C/+ equals to a contribution 
in U proportional to the first exponent, added by an additional contribution given by the first integral taken over 
a closed contour passing around the branch cut connecting the branch points. This correction, however, can not 
be represented in terms of elementary functions. Similarly, in the case p > b the first integrand possesses a pole 
ai z' = z + iO and a branch point at z' = i{b + p), while the second integrand possesses a unique singularity at 
z' — ~i{p — b). In this case, the branch cuts go from the branch points to imaginary infinity and do not intersect 
the real axis. Thus, the correction to the contribution in U proportional to the first exponent is given by two 
integrals over the banks of these cuts. Finally, at the ring p — b the second branch point approaches the real 
axis and becomes singular. Therefore, the nonsingular solution (|2ip is in fact a superposition of two singular 
counterpropagating beams. 

The spatial distributions of the electromagnetic energy density 

z^ = i-{[Re(E)f + [Re(H)f}, (26) 

in the three models of a monochromatic focused beam [the paraxial approximation defined by Eq. (jl3p . the 
ordinary complex source method defined by Eq. (I16p . and the modified complex source method defined by 
Eq. (PT|) ] are compared at the Fig. [T] in the tight focusing regime (6 ~ A). It is clear from the figure, that in this 
regime the field model corresponding to the ordinary complex source method occurs to be very close to paraxial 
approximation (the Narozhny-Fofanov model) everywhere but in the vicinity of the ring singularity, where it 
breaks down. At the same time, the solution obtained by the modified complex source method, being exact and 
nonsingular, differs from both of them. 
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Figure 1. The spatial distribution of the electromagnetic 
energy density in the monochromatic focused beam: the 
modified complex source method (curve 1), the complex 
source method (curve 2) and the paraxial approximation 
(the Narozhny-Fofanov model, curve 3). Panel (a): Dis- 
tribution in the focal plane (2 = <;!> = 0), panel (b): dis- 
tribution along the focal axis, and panel (c): distribution 
along the line parallel to the focal axis and crossing the 
ring singularity. The values of the parameters used are 
b = 1, go = 1, Lu = 5 {\ = 1.26). 



The same technique can be applied to construction of the pulse-like solutions as well. In order to do this, 
it is enough to introduce an envelop factor in q{t), e.g., of the Gaussian form e~* or of the Lorentzian 
form (1 -I- t^/T^)~^, where T is the duration parameter of the pulse. As it is clear from the consideration 
above, the result is an expression for the field components in terms of elementary functions, though generally a 
rather complicated one. Below, let us consider a special but perhaps the most interesting case of an ultrashort, 
subperiod focused pulse. An especially simple model of a pulse can be obtained if one drops the carrier frequency 
factor in q{t) totally. Furthermore, since for real values of t, z and p we have |Im(7^)| < h and hence —26 < 
Im(t — ib i: TV) < 0, in order to obtain a nonsingular solution for the whole range of the duration parameter, let 
us choose q{t) to be analytic in the lower half plane of complex t. Assuming q(t) — qo/{l + it/T) and applying 
(Pn| . we have 

The nonsingular solutions corresponding to single pulses propagating along (opposite to) z-axis can be extracted 
from Eq. (|27p by using ((24|) and are given by 

U±{v, t) = T , — =, (28) 

respectively. Thus, unlike the beam-like solution considered above, our solution ()27p represents a superposition 
of two nonsingular counterpropagating focused fields. A more general nonsingular solution of this kind can be 
constructed as an arbitrary superposition ci[/+(z, p, i) -f C2?7+(2:, p, — i) with complex coefficients. Note also, that 
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Figure 2. The distribution of the electromagnetic energy density in the xz-plane in subsequent stages of focusing of the 
subperiod pulse-like solutions corresponding to U (the top row), U+ (the middle row) and U- (the bottom row). The 
values of the parameters used are & = 1, T = 0.2. 



in the weak focusing limit 6 — > oo, T = const, and both expressions ([77|) and (|28p reduce to a plane wave moving 
along z with the profile U iq{t — z)/b. 

The resulting expressions for the vector potential and the electromagnetic fields are quite cumbersome and 
this short paper is not an appropriate place to adduce them. Rather, let us present a visualization of the 
electromagnetic energy distribution in the xz-plane, corresponding to Eqs. ([27]) and (|28|) . at different stages of 
focusing, see Fig. [21 It can be seen from the pictures, that even in the tight focusing regime U+ actually almost 
coincides to U (i.e., U- is negligible) except for the moments which are directly close to the moment of maximal 
collapse. Also, it can be recognized that unlike U and [/+, the backward propagating pulse J7_ possesses a 
coreless tubular structure. 



In conclusion of the section, let us present an independent derivation of the solution (|27|) on the basis of 
direct evaluation of the corresponding Fourier integral. According to the beginning of the section, the Fourier 
transform to start with is given by 



(29) 



Let us evaluate the Fourier integral U{r,t) = J c?a; d'^fc e*^'^ '"~'^*^C/ka;/(27r)* using the cylindrical coordinates in 
the fc-space. After performing the integration with respect to the azimuthal angle (which is trivial) and the 
integration with respect to (which is realized by taking away a (5-function), we are left with 



+ 00 Ul 

Uir,t) = 2tqoT ( duj e-'^^^'^'"-'^^ [ -^t^^ J„{k^p) cos[(z-z6)y^ 



(30) 



The inner integral over equals^"'^ sui{loTV) /TZ. Hence, the final integral over oj is elementary and we come to 
the formula (P7| . Note that in order to reproduce the Eq. ((25)) in the same way one should take the integral 
which is similar to the inner one in (|30p . but with the imaginary exponent instead of the cosine. However, we 
could not find it in the mathematical literature thus far. 



8 



5. CONCLUSIONS 

In this paper, we have expHcitly demonstrated that the "complex source method" ^^"^^ of generating the exact 

focused beam or pulse-like solutions of the wave or Maxwell equations should be modified in order to obtain the 
fields which are nonsingular in the focal plane, correspondingly to absence of the sources. Also, we have shown 
that this technique always generates solutions describing a pair of counterpropagating focused fields, rather than 
a single one. We suggest resolution to both disadvantages. In order to generate nonsingular solutions, one 
should exploit the difference between the retarded and advanced potentials instead of choosing a special branch 
of the retarded potential, as it was previously accepted in the literature. Besides, we propose the regular method 
which allows one to separate the resulting field onto the one propagating along the focal axis and the backward 
propagating one. However, even if one starts from the nonsingular solution, it is not guaranteed generally that 
the forward and backward propagating fields will remain nonsingular. Our final proposal is to start with an 
appropriate anzats for the vector potential in terms of scalar functions obeying the wave equation, and then to 
construct them by application of the complex source method. This allows one to control the polarization of the 
electromagnetic field beforehand. A reasonable example of such anzats with known polarization type (e- or h- 
circularly polarized wave) was proposed in Refs.'''^ 

As an application of the modified "complex source method", we obtained and visualized the solution which 
describes a counterpropagating focused monochromatic laser beams, and the solution describing a subperiod 
focused laser pulse. However, the technique under consideration is in no way restricted to the examples given 
and allows one to obtain a wide class of focused beam- or pulse-like fields. An advantage of this approach is that 
the field components can be expressed in terms of elementary functions, though generally in a rather complicated 
way. Nevertheless, it seems that such expressions are very suitable for numerical calculations of both the motion 
of the charged particles interacting with the laser field and the quantum effects induced by it, especially in the 
tight focusing regime. The key point for such applications of our expressions is absence of singularities in them. 
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